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Abstract 

Background: Salmonella enterica serovar Enteritidis is a zoonotic pathogen. Human infections are associated 
with contaminated eggs and egg products. In Japan, since 1989, the incidence of food-borne disease caused by 
5. Enteritidis has increased and a pandemic has occurred; however, little is known about changes that occurred 
before and after this pandemic event in the dominant lineage of isolates from food-producing animals. This 
study aimed to determine the 5. Enteritidis lineages in Japan over the last few decades by using multilocus 
variable-number tandem repeat analysis (MLVA). 

Findings: MLVA was used to analyse 79 5. Enteritidis isolates collected from chickens (n = 63), cattle (n = 12), pigs 
(n = 2), and goats (n = 2) during 1975-2009. The 5. Enteritidis isolates showed 14 different MLVA allele combinations, 
which were classified into two major clusters (A and C) and a minor cluster (B). All the 62 isolates in cluster A 
were isolated after 1988, whereas 13 of the 17 isolates belonging to cluster B and C were isolated before 1989. 

Conclusions: The MLVA results showed that cluster C was predominant before 1989, and isolates in cluster A 
disseminated since 1989 and replaced the previous dominant clone, suggesting that isolates of cluster A originated 
from imported 5. Enteritidis infection. 
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Findings 

Salmonella enterica serovar Enteritidis is a zoonotic patho- 
gen that can be transmitted to humans via many different 
reservoir hosts. Most outbreaks of S. Enteritidis infections 
have been associated with contaminated eggs and egg prod- 
ucts [1]. Until the 1980s, S. Typhimurium was the serovar 
most commonly isolated from humans worldwide; however, 
during the 1980s and 1990s, S. Enteritidis emerged as a 
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common cause of salmonellosis, first in European countries 
and then worldwide [2-4]. 

Similarly, in Japan, the serotype S. Enteritidis has been 
isolated most frequently, since 1989, from patients with 
food-borne illnesses [5]. Moreover, S. Enteritidis phage 
type PT8 was the most prevalent human isolate reported 
until 1988, until PT34 became the predominant phage type; 
PT1 and PT4 have replaced PT34 as the predominant 
phage types since 1992 [5]. However, little is known about 
changes that have occurred over the last few decades in the 
dominant lineage of isolates from food producing animals 
in Japan. 

Phage typing is a commonly used method for epidemio- 
logical surveillance of S. Enteritidis infection [6]; however, 
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Figure 1 Dendrogram of 79 Salmonella Enteritidis isolates based on MLVA profiles. The dendrogram is divided into three clusters: A, B, 
and C. The numbers in the columns SENTR-1 to SE-8 indicate repeat count. 
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it requires specialized reagents and laboratory equip- 
ment and does not always yield sufficient information 
for epidemiological purposes. Over the last decade, new 
techniques in molecular biology have been developed, 
such as pulsed-field gel electrophoresis (PFGE) [7,8], 2- 
enzyme ribotyping (Pstl-Sphl) [9], and multilocus variable- 
number tandem repeat analysis (MLVA) [10,11]. Of these 
methods, PFGE is now the gold standard for discriminating 
among strains at the DNA level [7]. However, based on 
amplification of a variable number of tandem repeat areas, 
MLVA is considered to have greater discriminatory power 
than PFGE and has been proposed as an alternative for 
genotyping of S. Enteritidis [10,11]. 

In this study, MLVA genotyping was performed to de- 
termine the lineages of S. Enteritidis before and after 
pandemics during the 1980s and 1990s. We have provided 
evidence that clonal dissemination and replacement have 
occurred among isolates from food-producing animals 
since 1989. 

A total of 79 isolates from chickens (63 isolates), cattle 
(12 isolates), pigs (2 isolates), and goats (2 isolates), which 



are epidemiologically unrelated, were collected in Japan 
during 1975 to 2009 in 24 prefectures in Japan (Figure 1), 
and these strains were used for MLVA. MLVA was con- 
ducted by amplification of 12 loci (SENTR-1, SENTR-2, 
SENTR-3, SENTR-4, SENTR-5, SENTR-6, SENTR-7, 
SE-3, SE-7, SE-4, SE-6, and SE-8) using previously designed 
primers [10,11]. PCR products were sequenced to confirm 
the number of tandem repeats. The number of tandem 
repeats at each locus was manually determined using 
Genetyx version 11.0 (Genetyx, Tokyo, Japan), and the 
motif numbers present in the tandem array were 
imported into BioNumerics version 6.0 (Applied Maths, 
Sint-Martens-Latem, Belgium). 

The 79 S. Enteritidis isolates included 14 different MLVA 
profiles (Figure 1). Nine distinct MLVA profiles were identi- 
fied from chicken isolates, 7 were identified from cattle, 2 
from goats, and 2 from pigs. A dendrogram was generated 
using the categorical coefficient and unweighted pair 
group method with arithmetic means by the BioNumerics 
software (version 6.0). MLVA profiles were classified into 
3 clusters (A, B, and C) delineated with a 67.8% similarity 
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Figure 2 Minimum spanning tree of 79 Salmonella Enteritidis isolates, (a) Each node represents a unique MLVA profile, indicated by the 
number in the middle of the circle. The nodal size is proportional to the number of isolates per MLVA profile. The colour within each node reflects the 
proportion of strains isolated from different periods per node. The length of the branches represents genetic distance (change in loci) between two 
neighbours. Bold short lines connect two MLVA types differing by a single MLVA locus. Thin longer lines connect double or triple variants, and dotted 
lines indicate connection between two types differing by four MLVA loci. Clonal complexes were created based on maximum neighbour distance 
changes at two loci, (b) This is the same minimum spanning tree as in Figure 1a. Wedges in circles indicate the proportion of isolates from respective 
sources with a particular MLVA profile. 
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cutoff value (Figure 1). Cluster A consisted of 7 profiles 
corresponding to 62 isolates, including 58 isolates from 
chickens, 2 isolates from goats, 1 isolate from cattle, and 1 
isolate from pig. Cluster B included 2 isolates from cattle, 
and cluster C, which was composed of 6 MLVA profiles 
representing 15 isolates, included 9 isolates from cattle, 5 
isolates from chickens, and 1 isolate from pigs. 

A minimum spanning tree (MST) was also created 
based on the categorical data sets (Figure 2). All the 
MLVA profiles in cluster A clustered into a single clonal 
complex in the MST, and all isolates belonging to cluster 
A were isolated after 1988 (Figure 2). With the exception 
of MLVA profile C6, profiles in cluster C were clustered 
into a single clonal complex, and this complex differed 
from the cluster A complex at 7 loci (Figure 2). The 
MLVA profile C3 included 5 isolates from cattle, 4 isolates 
from chickens, and 1 isolate from pigs (Figure 2). All 10 
isolates representing profile C3 were isolated before 1983 
(Figure 1). Of the isolates in cluster C, 4 isolates that 
exhibited MLVA profiles CI (from chicken) and C4, C5, 
and C6 (from cattle), were isolated in 1992, 1993, 1991, 
and 1997, respectively (Figure 1). These 4 isolates could 
have originated from isolates showing profile C3, which 
was the dominant lineage before 1983 (Figure 1). A profile 
of cluster B, comprising 2 isolates collected from cattle in 
1982, exhibited changes in the neighbour distance at 3 
and 4 loci from clusters A and C, respectively (Figure 2). 
Therefore, the lineages of these 2 isolates from cattle seem 
to be different from those of clusters A and C, which 
included isolates from chickens. A proportion of the 
cluster A isolates isolated after 1988 was statistically 
different from that of isolates before 1989 (P < 0.01, by 
Fishers exact test). Taken together, these results indi- 
cate that MLVA cluster C was predominant before 1989 
and that clonal dissemination of lineage MLVA cluster 
A has occurred in food-producing animals in Japan 
since 1989. Porwollik et al. [12] differentiated between 
two major phage type lineages, including PT4 and PT8, 
by whole genome microarray analysis. These lineages 
might correspond to MLVA clusters A and C, respectively. 
Since PT8 was the most prevalent phage type until 1989 
in Japan [5], isolates of this phage type might belong 
to MLVA cluster C. S. Enteritidis phage type PT4 was 
detected in isolates from chickens imported from England 
in 1989 [8], and a variety of phage types, including PT4, 
was observed in isolates from chickens after 1989 [8,13]. 
Although we have not yet determined the phage type 
and have not differentiated between imported and 
native isolates, no isolate classified into cluster A was 
detected before 1989, suggesting that isolates in cluster 
A may have originated from imported S. Enteritidis in- 
fections; consequently, clonal expansion of S. Enteriti- 
dis might have occurred in chickens and also affected 
humans. 



Since MLVA profiles can be easily stored in a database, 
routine and long-term epidemiological surveillance with 
such a method may enable early recognition of potentially 
epidemic Salmonella clones. 
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